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, The functional integral representation for fermionic observables on the lat- 

I tice is studied. In particular, Grassmannian representations of the scalar J^"^-* 

and pseudoscalar J^^^ currents and pseudoscalar correlator are derived. It is 
^ ' also discussed the connection between the fermionic Fock space and boundary 
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conditions along the time direction. 



1 Introduction 



^ I The 'functional' integral approach to the quantization with lattice regularization in the 

euclidean space is very convenient for numerical calculations and gives a possibility for 
the non-perturbative study of the gauge theories [0. However, a consistent quantiza- 
tion scheme needs a canonical Hamiltonian (or transfer matrix) approach to answer 
questions that do not obviously hold in the functional integral formulation. 

The canonical quantization formalism is necessary to define boundary conditions 
for and ip,ipx in the integrals in eq.'s (|1.5[),(P^), in particular along the 

forth (imaginary time) direction. It is the canonical quantization formalism which 
establishes the connection between correlators of currents (e.g., pseudoscalar current 
J^^^) and masses of corresponding particles. 

The standard Wilson action Sw with SU{Nc) gauge group is ^ 

Sw{U; ?) = Sg{U) + Sp{U; ^) , (1.1) 
where pure gauge part Sg and fermionic part Sp are given by 
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f=l x,y s,s'=l 
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In equations above U^^ G SU{Nc) are gauge fields, ipxytpx fermionic Grassmann 
variables and g = go{a) is a bare coupling, a being a lattice spacing. Wilson's fermionic 
matrix M.{U) (with Nf = 1) is given by 



11,11 



:i.3) 



where k is a hopping parameter and -Pj^-* = (1 ± 7^)/2 . The action Sw is invariant 
with respect to the local transformations 



Uxn — > Uxi, = ^xUx^l^l+^ ; tpx ^xipx ; ^x '^x^l , (1-4) 
with fix £ SU (Nc) . The partition function Zw is 

Zw = j[dU][di)(kp] e-5^(^''^'^) , (1.5) 

where [dU] denote product rix^t dUx^ , etc.. The average of any functional 0{U ; V'; iIj) 
is chosen to be 

(C) = ^ f[dU][d^di^] 0{U]i)]i^) ■ e-^>^(^;^'^) . (1.6) 
Zw J 

Given some boundary conditions, the average {O) is (mathematically) well defined 
and can be calculated numerically. 

In the canonical quantization approach, the connection between the transfer matrix 
V and the Hamiltonian H is given by 

V = e-"^ , (1.7) 
The corresponding partition function Zh is 

= Tr fe-^^) =Trfr^4) , (1.8) 

V J colorless ^ ^ colorless 

where A^4 is the lattice size along the forth direction, T = l/aN^ = I/L4 and the trace 
is defined on some colorless space of states. The consistency between two definitions 
of the partition function given in eq. ( |1.5| ) and eq. (pTS]), i.e. Zw = Zh = Z , defines 
the transfer matrix V 0, |], Q . The average of any field operator O is 

id) = ^Tr (dV^A . (1.9) 

^ ' Z ^ /colorless ^ ' 
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Let be eigenstates of the transfer matrix V with eigenvalues 

V\^k) = Afcl^fc) ; Xk = 6-'^^'= , (1.10) 
and Ek are eigenvalues of the Hamiltonian. Then 

(a)~5:e-^^^-(M/fciaiM/fc). (1.11) 

fc>0 

Similar expressions can be written for the correlators of currents, e.g. pseudoscalar 
current J^^\ etc.. 

The question of interest is the connection between operators O defined as nor- 
mal products of creation and annihilation operators and corresponding functionals 
0{U;ip;ip) in the functional integral approach. Another problem of interest is the 
choice of the boundary conditions along the imaginary time direction. 

This paper is dedicated to the connection between the operator (canonical quan- 
tization) formalism ^] and functional integral approach The transfer matrix 
formalism is briefly reviewed in the second section. In the third section the statistical 
averages of the scalar J^^^ and pseudoscalar J^^^ currents are calculated as well as 
pseudoscalar correlator. Boundary conditions are discussed in the forth section. The 
last section is reserved for conclusions. 



2 Transfer matrix formalism 

Let us give an outline of the transfer matrix formalism following the paper ||^. The 
main modification is connected with the introduction of the projection operator Pq (see 
below) which is necessary to take into account the Gauss law ^. 

Let c|(cj) and dl{di) be creation(annihilation) operators of fermions and antifermions, 
respectively, that satisfy canonical anticommutation relations : 



dj 



(2.1) 



and other anticommutation relations are equal to zero. Indices i,j are composite : 
i = {x; a; s) where x is a three dimensional coordinate, a = 1, . . . , N^. is a color index 
and s = 1, 2 is a spin index. Therefore, i,j = 1; . . N where = 2A''cV3 and V3 is 
a threedimensional volume. 

Let {Ug-k} and {Ug.f,} be two configurations of the gauge fields defined on the 
spacelike links Ig = {x; k) . The transfer matrix V is an integral operator with respect 
to the gauge fields. Its kernal V{U; U') is an operator in the fermion Hilbert space 
and has the following form P| 



ViU;U') = fl{U)VG{U-U')fF{U')- (2.2) 



^In paper the gauge U^-a — 1 has been chosen which is not possible on the finite torus. 
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fpiU) = Co-exp[d^Q{U)c}exp[-c^R{U)c-d^R^{U)d} ; 

fliU) = Co-exp[-c^R{U)c-d^R^{U)d}exp[c^Q{U)d^^} , 

where Vg{U]U') = {U\Vg\U') corresponds to a pure gauge part 0. Hermitian 
matrices Q, R are given by 
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and 
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Qxy ~ 2 ^ [^y,x+kUxk - Sy,x-kUl_i, ,, ■ CJfc . (2.4) 

fc=l 

Co = C'o(f^; U') is a constant to be defined later and ak are Pauli matrices. 
It is convenient to define Grassmannian coherent states 



|r/C) = exp + 4Cx) } |0) ; {vC\ = (0| exp + Qd,) | , (2, 



5) 



where rjg, . . . ,(x some Grassmannian variables (spin and color indices suppressed). 
It is easy to see that 

cM) = vM) ; (^Cl4 = (^Cl% , (2.6) 

and 

(VCWC) = 6^^^'+^^' , i = J[dr]dv][dCdC] e-^^-~<^\7^0{vC\ ■ (2.7) 

To prove the equivalence between the transfer matrix and functional integral ap- 
proaches one can start with 

Zh = Tr (t/^*) = Tr (t/^^Po) , (2.8) 

\ / colorless \ / 

where Pq is the projection operator on the colorless state 

Po = Jm 7^(A) ; [dA] =lldA{x) , (2.9) 

X 

and 7^(A) is the gauge transformation operator. In particular, 

7^(A)x^7^t(A) = Alxx ; n{A)xln\A) = xIa, , (2.10) 



where 



x.= ( ; xt = (4 4). (2.11) 

To obtain the functional integral representation of Zh one needs to represent as 
a product : V'^'^ = V -V ■ . . . -V and insert unit operators 1 defined in eq. ( |2.7| ). For 
every time slice one can introduce new Grassmannian variables ip£{x4),'ipg{x4) : 

M^^) = E B.;/^ f ^1 ) ; Mx.) = E(% - Cj)Bj/^ . (2.12) 

Therefore, [drjdrj\[dC,dQ = ■ [dipdip] . The constant Cq is chosen to cancel the 
Jacobian J of this transformation. Assuming that the fermionic Fock space spanned 
by all possible fermionic states, one obtains 

r ^ — — 

Zh= l[dA^l[Y[dUsk{x4)Y[(^^xdi^. exp{-Sw{U-,^-,^)}, (2.13) 

where Sw is the Wilson action with f/f4(x4) = 1 and boundary conditions 
tp^{xi = N^a) = -ip§(x^ = 0) ; 

f/xfc(x4 = N^a) = f/4(x4 = 0) . 
The final step is a change of variables 

i^xi^i = N^a) — > Agipg{x4 = N^a) ; V'x(2^4 = N^a) — ^ i'xi^i = ^4*^)^1 > (2.15) 
and similar for the gauge fields U. One obtains 

r ^ _ _ 

Zh = YldllMix^ = N^a - a)YlYldU:gk{x4:)[dtljd'ilj] exp \^-Sw{U;iIj;iIj)^ , (2.16) 

X X k = l 

where 

U^4 = l at < 2:4 < A^4a - 2a ; Um{xa = N^a - a) = , (2.17) 
and boundary conditions 

i'xixi = N^a) = —ip^^x^ = 0) ; i'xi^i = N^a) = —i^^ix^ = 0) ; 

(2.18) 

U£k{x4 = N^a) = U£k{xi = 0) . 



ip^ixi = N^a) = -i/j^{xi = 0) ; 

(2.14) 
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Evidently, Zh coincides with Zw in the temporal gauge defined in eq. ( |2.17| ). This 
proves the equivalence of two approaches. 
Two comments are in order. 

i) In the functional integral formalism Grassmannian variables ?], rj, C, C correspond 
to the operators c,c\d,S. However, new Grassmannian variables ip^ip are con- 
nected with TjjTjXX in ^ rather nontrivial way given in eq. (|2.12| ). This ob- 



servation will appear to be important for the calculation of fermionic matrix 
elements. 



ii) Antiperiodic boundary conditions in eq. ( p.l8|) stem from the choice of the fermionic 



Fock space which is, in fact, a physical assumption. Another choice of the 
fermionic Fock space gives another boundary conditions for fermionic variables 
ijj,'ip along the time direction. In details this question will be discussed later. 



3 Fermionic currents 

3.1 Pseudoscalar current 

Pseudoscalar current is given by 

JP = -ixh^l^Xx- = iicidi'^ - d^c^) , (3.1) 
where 74, 75 are euclidian 7-matrices 

Ik = ( T ) ' = ( -1 ) ' ■ (2-2) 

Evidently, J^'^ = J? . Let us derive the functional integral representation of the 
statistical average 



= ^Tr (J(^) V>^^Po) , (3.3) 



where 



J(P) ^ ^ jiP) . (3.4) 



As well as in the previous section, we assume that the fermionic Fock space spanned by 
all possible fermionic states. To derive the desired expression one should proceed in a 
same way as outlined in the second section and use the properties of the Grassmannian 
coherent states in eq. (|2.6|) . Taking into account 



(%G - Cjw) (3^4) = ^27(3^4) Bi/,^(x4)75BVi(x4)V^^-(x4) , (3.5) 

yz 

one obtains 



6 



= _^ j[dU][d^^] (^B75V^)(0) ■exp{-^H/(f/;V^;V^)} , (3.6) 

where 

(V'B75'0)(O) = X!^x(a^4)B^^(a;4)75^y(a;4)|^ , (3.7) 



and boundary conditions given in eq. ( p.l8| ). 

For the zero-momentum pseudoscalar correlator 



r(^)(r) = ixr (K^^-V(^)rV(^)Po) , (3.8) 



one arrives at 



r(^)(r) = -i f[dU][dijij] [^pBJ,i:){T)(^pBJ,iJ){0) ■ exp{-Sw{U;ij-ij)} . (3.9) 



One can see that Grassmannian current J^^'^ has no naive expression J^i^e — 
J2x4'xl5'^x ■ Instead, it has rather comphcated expressions and depends on the gauge 



field Uxfj,- However, in the continuum limit a 



_ ^ _ _ 

V'x75^x (i^xUxklbi^x+k + V^^f/i.fc^fcTS^/'x-fc) (2:4) 

k=l 



^ (l-6«:)X:(^.^75^.^)(a:4) ~ 42.e • (3-10) 

X 

Therefore, at nonzero spacing a there is a nonperturbative renormalization of the 
matrix element and correlator. 



3.2 Generalized partition function and chemical potential 

Let us derive the functional integral representation for the generalized partition func- 
tion Z{\q, \q) defined as 

Z(A„ A^) = Tr (t/^^ ■ exp{A,iV, + \%}Po) , (3.11) 

where 

% = Y. ; % = 11 4dx . (3.12) 

X X 

Choosing the standard fermionic Fock space and repeating calculations outlined in the 
second section one obtains 

Z(A„ Xg) = J [dU] [d^dtfj] exp {-Sw + SSp} (3.13) 



with 



5Sp 



eA,+A, _ i)J2M0)Pt^C,gmm , (3.14) 



xy 



and 



1 ^ 



t 



7A: 



(3.15) 



fc=i 



The partition function 2^(/i) with nonzero chemical potential fi is given by 

Z(/.)=Tr (exp{-i(]f-/x]V)}Po) =Tr (^t/^* . exp{|]v}Po) , (3.16) 



where N = Ng — Nq . Choosing 



one obtains 



A - A - ^ 



(3.17) 



Making the change of variables 



(3.18) 



e-^4MV's(x4) X4 7^ 



and '^^{x^ 



one obtains the modified fermionic matrix A4{U; fi) 



e^'4^^^(x4) Xi^O 
e^'*^^^(x4) X4 = 



(3.19) 



t 

x—k,k 



k=l 



-2k 



^ "s/;a;+4-' 4 ^ xA ^ c Oy;x-4-' 4 '^a;-4-4 



(3.20) 



Evidently, Ai{U; fi) coinsides with the fermionic matrix for the nonzero chemical 
potential proposed in 
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3.3 Scalar 'condensate' 

Let us obtain the functional integral representation of the statistical average {J^^^) 
of the scalar current 



In what follows this average will be refered to as scalar condensate. Evidently, 

1 d 



(J(^)) = iTr (J(^) t/^^Po) = ^^^(A) 



A=0 



where 



Z(A) = Tr (t/^^ ■ exp{A(iV, + N^)}Po) 
From eq.'s ( p.l3| ),( prT4D one derives at Xg = Xg = X 

Z{X) = J[dU][(hpd^] exp{-Sw + SSf} ; 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



6Sf 



2k 



-(e^^ - l)^(0)Pi"^C(0)^(0) 
where the notations are used 



ij{a)Pt^Ul{0)^{0) = EV^a-(^4 = a)Pt^Uyx4 = O)^^^ = 0) 

X 

etc.. Finally, one obtains 



(3.25) 



(J(^)) = 1 [dU][di)di)] e-^'^^^'^'^) 



Z 



■2k 



\p{a)Pt^Ul{0)iP{0) + ^(0)Pi~^?74(0)^(a) - 2^(0)Pi ^C(0)^(0) 



Jile = (^V^)(0) = = 0)V^-(^4 = 0) . 



(3.26) 



Evidently, the Grassmannian current 

J(^) = 2K\^{a)Pt'^Ul{0)ij{0) + V^(0)Pi"^f/4(0)7/'(a) - 2^(0)Pi"^C(0)V^(0)] , (3.27) 
does not coincide with 



(3.28) 



In the continuum limit a ^ one obtains 

j(s)_2«:(^^)(0)~41t • (3.29) 

As well as in the case of the pseudoscalar current, there is a nonperturbative renormal- 
ization of the scalar condensate at nonzero spacing a . 



4 Fermionic Fock space and boundary conditions 

One or another choice of the fermionic Fock space depends on the model (physical) 
assumptions. For example, QCD vacuum is expected to have an equal number of quarks 
and antiquarks and it puts a restriction on the choice of the states. The important 
observation is that boundary conditions for Grassmannian variables ■?/', ifj along the 
imaginary time direction depend on this choice. 

Till now it has been assumed that the fermionic Fock space is spanned by all possible 
fermionic states, i.e. states 



i{M;w.}) = n(4rn(4no), (4.1) 

where ni,mj = 0,1 and = 2NcV3 is a maximal value of parameters i,j . This 
choice of the Fock space results in the antiperiodic boundary conditions for fermionic 



Grassmannian variables in eq. ( |2.18| ). Indeed, the trace of any fermionic operator O 
can be written as 



Tr^O = J[dr]dr]][dCdC][didr]'][dC'dC] exp{-r]r]- a- r]'r]' -('('} 

■{vC\OWO-TT{\r^0{v'C\) ■ (4.2) 
Choosing Fock space as in eq. (|4.1|) one obtains 

Tr {ivOWCl) = exp i - j:{rf,V^ + ZCi)] , (4.3) 



4 = 1 



and 

T^r^O = J [df]dr]] [dCdQ expj-rjr^ - CC} " (^C|0| - r/; -() . (4.4) 

Minuses in \ — rj; —() presume antiperiodic boundary conditions for ipjip along the 
imaginary time direction. 

However, the choice of the fermionic Fock space made in eq. ([4.1|) is not a unique 
one. As an example, let us consider the zero temperatures limit T 0. In this case the 
main contribution to the partition function Z is expected to come from the vacuum 
eigenstate |vac) : 
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Z = (vac|t/^^|vac) ; T . (4.5) 

The corresponding Fock space is supposed to have the equal number of quarks and 
antiquarks. Therefore, it looks reasonable to choose the fermionic Fock space spanned 
by the vectors 

N N 

|{n.};{m,}) = n(4rn(4ri0) ^Ith ^n, = ^m, . (4.6) 

i=l j = l i i 

Using the following representation for the delta-symbol 

= ^ E exp I '^kY,{n, - m,) I , (4.7) 



fc=0 1 ] 



one obtains 



1 _ 

Tr [\riO{vC\) = ]7 E exp {-qiJi'v - qK' C) , (4i 



fc=0 



where 



gfc = exp|^| ; A: = 0;l;...;iV-l. (4.9) 



Therefore, 



1 ^"-^ r - - 

TrxO = ^ E J[dvdv][dCdQ expj-rjr^ - CC} • {vC\0\ - qkV, -qIO ■ (4.10) 

k=0 

In the infinite volume limit V3 ^ oo one arrives at 

Tr = i- j[dr)dri][dCdC] exp{-rjr/ - CC} • (r/C|0| - e^^r/; -e'^^C) • (4-11) 

Finally, fermionic boundary conditions for ip,ip are 

^^[xi = N^a) = -e''^ip^{xi = 0) ; V'x(x4 = N^a) = -e'"^i^g{x4 = 0) (4.12) 

with integration over if between and 2tt. On the finite lattice the temperature is 
nonzero and higher exitation states can give the (artifact) contributions. Therefore, 



one may expect that at A^4 < cxd the boundary conditions given in eq. ( |4.12| ) could be 
a better choice for the zero temperature calculations (e.g., for the hadron spectroscopy 
study) as compared to the choice of antiperiodic boundary conditions. 

Another interesting case is the finite temperature transition(s) in the Universe. The 
baryon asymmetry AB of the Universe is small and it is expected to be zero in QCD. 
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Therefore, also in this case the choice of the boundary conditions given in eq. ( [4.12| ) 
looks preferable. It is worthwhile to note that for Polyakov loop V one obtains 

(P)=0, (4.13) 

and (1^1) is expected to be a good order parameter as in quenched QCD. 

Of course, different physical problems need different assumptions about the struc- 
ture of the fermionic Fock space. For example, finite temperature phase transition with 
nonzero baryon density (Ai? 7^ 0) presume another choice of the fermionic Fock space. 
One can expect that the phase transition temperature Tc depends on the choice of this 
Fock space and, therefore, on the choice of the boundary conditions along the forth 
direction. 



5 Conclusions 

The first goal of this work is the study of the connection between the fermionic currents 
J in the canonical quantization approach and corresponding currents J in the func- 
tional integral (Wilson's) scheme. In the canonical quantization approach the fermionic 
currents J (e.g. pseudoscalar current, etc.) are defined as normal products of oper- 
ators Xx, Xgy while in the functional integral formalism the fermionic currents J are 
constructed out of Grassmannian variables ipxi'^'x- 

As an example, two operator currents J^^^ = :xi74Xx- and J'^^^ = :xl;l4l5Xx'- 

have been considered. It appears that corresponding Grassmannian currents J have 

(s) — (p) — 

no naive expressions JnaLe = i^xi^x and JnJve = '^xl^'^x ■ Instead, Grassmannian 

currents have rather complicated expressions and depend on the gauge field Ux^j.- How- 
ever, Grassmannian currents J^^'^^ become ~ J^ai^l in the continuum limit a — > . 
Therefore, at nonzero spacing there is a nonperturbative renormalization of the matrix 
elements and correlators. 

Another scope of this paper is the study of the boundary conditions along the 
imaginary time direction. The important observation is that boundary conditions for 
ipjip along the imaginary time direction depend on the choice of the fermionic Fock 
space, and the choice of the fermionic Fock space depends on the model (physical) 
assumptions. For example, in the zero temperatures limit the main contribution to the 
partition function Z is expected to come from the vacuum eigenstate |vac) which is 
supposed to have an equal number of quarks and antiquarks (if any). In this case it 
is reasonable to choose the fermionic Fock space as shown in eq. ( |4.6| ) which results 
in the boundary conditions for ipjip given in eq. ( [4.12|) . One may expect that on 



the finite lattice with A^4 < 00 the choice in eq. ( [4. 121) can be better one for the zero 
temperature calculations as compared to the case of antiperiodic boundary conditions. 
These boundary conditions look also reasonable for the study of the finite temperature 
transition on the cosmological scale with zero baryon density. 

In this paper the derivation of Grassmannian currents has been performed in the 
theory with Wilson fermions. It could be interesting to make a comparison with another 
versions of lattice QCD, e.g., with staggered fermions 0. 
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